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1 Introduction 


The goal of this grant was to study the properties of specific wide-band-gap 
semiconductor alloy systems of interest to the Electronics Division at ONR. 
More specifically this ‘grant has focused on alloys related to the group-III 
nitrides. This class of materials is of great relevance to a variety of electronic 
and opto-electronic devices, notably blue LEDs (light emitting diodes) and 
LDs (laser diodes). The alloy properties play a crucial role in band structure 
engineering of the materials used in such devices. 

The general questions of interest on semiconductor alloy systems which 
were investigated in this grant are: 


e miscibility and phase diagrams 

e band-gap bowing and direct vs. indirect transitions 
e ordering effects 

e crystal structure effects (e.g. wurtzite vs. zincblende) 


Specific systems under study were Al,B,_,N, GaN,Pj_2, IngGa,_,N, In, Alh_;N, 
and Al,Ga,_,N. We also investigate a compound LiGaQ, because we con- 
jectured that it could be viewed as a model system for cation ordering in 
wurtzite III-N alloy systems. It is in addition of interest as a closely lattice 
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matched substrate for GaN growth. The references below refer to the list of 
publications in section 0. 


2 Approach 


The approach we followed was based on first-principles (i.e. local density 
functional) electronic structure calculations using the full-potential linear 
muffin-tin orbital method. This method provides the basic properties for 
well-defined crystal structures. To describe the disorder effects in the alloy 
systems, the methodology followed was to describe them in terms of a statis- 
tical average of the corresponding properties in particular crystal structures. 
Since these represent the various local environments, this technique is known 
as the cluster expansion method. Details are provided in our published pa- 
pers. 

A crucial question is how to deal with the bond length relaxations in those 
systems which contain atoms of different size. A good deal of our investiga- 
tions was concerned with this question. Originally, we had hoped to be able 
to handle this with a semi-empirical approach such as the Keating model. 
This however proved to be insufficient and it was found that it could lead to 
inaccurate predictions such as deviations from Vegard’s law [1], in particular 
in the way we had implemented it at first, in which the crystal was first re- 
laxed vs. volume with ideal positions and only subsequently the atoms were 
allowed to relax at fixed volume. Subsequently, we spent a considerable effort 
to fully relax each of the structures under study to their minimum energy [2]. 
Finally, we realized that in a true random alloy, there is a certain amount 
of residual strain which has to be incorporated. We recently formulated a 
simple approach based on the actual bond length distribution as a function 
of concentration to handle this problem. This approach is similar to one 
that we previously had successfully applied to metalic alloy systems such as 
NiPt.’ 

1¢, Amador, W. R. L. Lambrecht, B. Segall, and M. van Schilfgaarde, Strain effects 
on the phase diagram of Ni-Pt alloys, Phys. Rev. B 47, 15276 (1993). 


3 Technical Results 


The most important results obtained during this grant are: 


1. The electronic band structures and relaxed total energies were obtained 
for a set of zincblende based crystal structures for the In-Ga, In-Al and 
Al-Ga nitride alloy systems. ‘These structures include the purely bi- 
nary end compounds, the 50 % compounds in three different crystal 
structures: [1g (or (001)-ordering), £1, (or (111)-ordering) and chal- 
copyrite, and the 25 % structures in the L1> structure. 


2. From these, the energies of formations as a function of concentration 
were determined for the disordered completely random alloys and com- 
pared to those of ordered structures. Generally, the disordered struc- 
tures are found to have lower energy of formation than the so-called 
Connoly Williams structures of the same concentration. However, the 
1, structure has higher energy and the chalcopyrite structure is a low 
energy structure. These findings agree with other alloy systems such 
as e.g. InGaAs. 


3. From the band structures as a function of concentration, we can pre- 
dict the direct-to-indirect transition in Al-containing alloys systems. 
In the Al-rich limit the gap switches from being direct at I to indirect 
at X. We also obtained the bowing coefficients. Finally, we obtained 
an interpolation formula for the band-gap bowing in the fully (cation) 
ternary system (or in more conventional terminology the fully quater- 
nary system) Al,Ga,In,N with z+y+2z=1. 


These results were published in [2]. 


4. In the last stages of the proposal, we developed an “effective volume” 
approach to predict the miscibility phase diagrams. ‘In this approach, 
we assume that the net effect of long-ranged elastic interactions leads 
to a distribution of bond lengths characterized by a certain amount 
of residual strain. We assume more precisely that the various nearest 
neighbor bond lengths vary linearly with concentration and can be de- 
termined by examining the dilute limits. For example, to determine the 
bond lenth of In-N in In-Ga alloys we examine the In-impurity in GaN 
and In-N bondlength in pure InN. The latter can simply be obtained 


from a Keating model as described by Chen et al.2 We use Keating 
model parameters obtained from our studies of the elastic constants.? 


Next, we assume that in an alloy at concentration x, the nearest neigh- 
bor tetrahedron clusters of each type (4 In, 3 In and 1 Ga, etc.) have the 
corresponding bond lengths rather than the bond length which would 
minimize this structure when found in isolation in a pure crystal. This 
assumption reflects the frustration in finding the optimum relaxation 
in a disordered alloy. ‘This leads effectively to a larger energy of forma- 
tion than assuming perfect relaxation and thus quantifies the residual 
strain effect on the alloy energetics. From the complete set of results 
of the various crystal structures as a function of volume, we extract 
cluster energies (x) as function of bond length and hence for each 
concentration we can simply describe the energetics as a weigted av- 
erage of E,,(x) applicable at that particular concentration. Combined 
with the regular solution model we can then determine the spinodal 
and binodal curves of the miscibility gap. Our miscibility gap temper- 
atures for InGaN are quite high, consistent with the observed tendency 
towards phase separation in this system. More quantitative details will 
be provided in a final paper, which is in preparation. 


5. Some cation ordered wurtzite based structures were also calculated and 
their internal parameters relaxed. Their complete analysis in compari- 
son with the zinblende results remains to be carried out. We find the en- 
ergetics for wurtzite to be similar to those in corresponding zincblende 
alloys. In other words, the local chemistry (number of In or Ga neigh- 
bors to a N) is more important than the stacking (cubic or hexagonal). 
Some of these structures were inspired by the structure of LiGaQOnz, 
discussed below. 


*A.-B. Chen and A. Sher, Semiconductor pseudobinary alloys: bond-length relaxation 
and mixing enthalpies, Phys. Rev. B 32, 3695 (1985). 

3K. Kim, W. R. L. Lambrecht and B. Segall, “Elastic Constants and Related Properties 
of the Group-III Nitrides,” in Gallium Nitride and Related Materials, p. 399 edited by 
R. D. Dupuis, F. A. Ponce, S. Nakamura, and J. A. Edmond, Mater. Res. Society 
Symp. Proceedings Vol. 395, (MRS, Pittsburgh 1996). and “Elastic constants and related 
properties of tetrahedrally bonded BN, AIN, GaN, and InN,” Phys. Rev. B 53, 16310-26, 
(1996); and Erratum Phys. Rev. B 56, 7018 (1997). 


6. We studied the band structure and total energy properties of LiGaQ, 
in various structures.[3] We find the chalcopyrite structure based on 
zincblende to be close in energy and properties to the actual crystal 
structure, while alternative orderings in the wurtztie based lattice have 
very different energy and band gaps. We also determined the band off- 
set of this material to GaN using the dielectric midgap level approach. 


4 Impact on other research 


The grant was complemented by an experimental AASERT award which 
had as original goal to grow Al,B,_,N alloys. In the later stages of this 
AASERT grant, a novel method was developed to grow GaN from liquid 
Ga and InN from liquid In with a plasma activated nitrogen source. As 
part of these activities, we studied the phonon spectra of InN in order to 
support the experimental characterization by Raman spectroscopy [4]. This 
work also formed the basis for our renewal proposal which will study of the 
present grant which was meanwhile awarded and will among others study 
the growth process of GaN from liquid Ga. 

The support trough this grant also allowed the PI to devote time to 
writing of two review articles on the band structures of the nitrides [5,6]. 
He also collaborated with a group of North Carolina State University, on a 
paper on SiC/AIN band offsets which was supported in part by ONR [7]. 
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BAND-GAP BOWING IN Ga,_,JIn,N ALLOYS 


WALTER R. L. LAMBRECHT 
Department of Physics, Case Western Reserve University, Cleveland. OH 44106-7079. U.S.A. 


Abstract—The band-gap bowing in Gai_,In.N is studied by means of linear muffin-tin orbital band 
structure calculations within the local density functional approach. A ‘‘virtual crystal”’-like approximation 
is compared with the Connolly—Williams cluster expansion approach with and without volume and 
bond-length relaxations. A progressively stronger bowing is found as the physical model's complexity is 
increased. The final bowing parameter obtained is 2.6 eV, which is significantly larger than the value 
reported in the literature. © 1997 Elsevier Science Ltd 


1. INTRODUCTION 


The band-gap bowing in Ga,_,In,N alloys is of 
significant interest for their use as active light- 
emitting layers in opto-electronic devices because the 
band-gap emission in this alloy system spans the 
violet to orange regions of visible light. Only limited 
experimental information[l] is available on the 
band-gap bowing because of the difficulties in 
obtaining homogeneous and high quality Ga, _,.In,N 
of high In content. From the theoretical point 
of view, the significant lattice mismatch between InN 
and GaN (about 10%) makes this system more 
challenging than the Al,Ga,_,.N studied in earlier 
work{2]. 

In this paper, we present a study of the band-gap 
bowing with a progressively more complete physical 
model. We progress from a virtual-crystal type 
approximation to a cluster expansion statistical 
description within ideal structures, to a model in 
which volume relaxation and finally bond-length 
relaxations are taken into account. Ordering effects 
on the band gap are also briefly touched upon. We 
find a progressively stronger bowing as the complex- 
ity of the model 1s increased. After briefly describing 
the underlying computational method, we describe 
the results and end with a comparison of our results 
with literature data. 


2. COMPUTATIONAL METHOD 


The first principles band structure calculations 
reported here were carried out using the linear 
muffin-tin orbital (LMTO) method{(3] in the atomic 
sphere approximation (ASA). The potential in this 
method is determined self-consistently within the 
local density approximation (LDA) to density 
functional theory(4]. As is well known, the LDA 
Kohn-Sham eigenvalues, which are not quasi-par- 
ticle eigenvalues, underestimate the band gaps—by 
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1.0 eV for GaN and 1.7 eV in InN[5,6]. However, the 
required self-energy correction is expected to be fairly 
independent of structure{7,8] and can thus simply be 
linearly interpolated for the alloys. The calculations 
reported here are for the zincblende structure. We 
note that the valence-band maxima and conduction- 
band minima in both GaN and InN and in both 
zincblende and wurtzite occur at the I point. Since 
the corresponding states are very similar in nature in 
both structures, we expect the gap bowing in 
zincblende to be nearly the same as that in wurtzite. 
The only difference is that the absolute values of the 
gap are slightly smaller in the former—by 0.4 eV in 
both InN and GaN. The details of the methodology 
for alloy modelling are described below. 


3. RESULTS 


The first model we applied is a virtual crystal type 
approximation. Within the LMTO-ASA method, the 
Hamiltonian is described by a set of potential- 
independent structure constants Spe. where R 
denotes atomic sites and L=(l/,m) the spherical 
harmonic component, and a set of potential 
parameters—C,, Aa, yu and p,-—describing the 
scattering at each type of atomic sphere a and for 
each angular momentum channel / = (s,p,d). The Cy 
define essentially the center of the corresponding 
unhybridized bands which can be thought of as 
atomic energy levels, and the A, define the 
corresponding widths ofthe bands and can be thought 
of as parametrizing the hopping integrals as 
J AwSeiee/Arr. We now define virtual crystal 
averaged potential parameters by interpolating the 
corresponding parameters for InN and GaN, i.e. 


Cu(x) = xCP + (1 — x)CH, (1) 
and 


Aux) = (Az)(As)'—*. (2) 
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The remaining parameters define further details in the 
dispersion of the bands which are less atom specific. 
To be unbiased. we average them linearly as in 
eqn (1). Within this approach, each cation is defined 
as an average cation with mixed In and Ga-like 
scattering properties and each nitrogen, N, is defined 
as an effective N with average scattering properties of 
those in InN and GaN. (We also average the 
parameters of the empty spheres that are introduced 
in the interstitial region as is necessary for open 
structures like zincblende.) This LMTO virtual 
crystal approximation (VCA) is in fact closer in spirit 
to the average f-matrix approximation because it is 
the scattering properties that are averaged rather than 
the potential. The results of this approach are shown 
as the squares interpolated by the full line in Fig. 1 
and show a fairly smail band-gap bowing. For 
convenience the linear interpolation between GaN 
and InN is indicated by the long dashed line. 
Next, we incorporate the effect that in a real alloy, 
the atoms are actually either In or Ga rather than an 
average atom and thus different local neighborhoods 
exist around each N atom. This can be done by means 
of a cluster expansion approach(9]. Considering the 
tetrahedral nearest neighbor environment of each N 
atom, its surrounding cations can either be 4 Ga, 
or 3Ga and 1 In, or 2Ga and 2In, or | Ga and 
3In, or 4In. In the alloy of a given average In 
concentration x, each of these cluster configurations 
which can be specified by the number of In atoms 
n=0,..., 4 occur with certain probabilities P,(x,7) 
which generally depend on temperature 7. In a 
general cluster expansion technique, not only 
tetrahedral nearest neighbor clusters but longer-range 
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Fig. 1. Band gap of zincblende based Ga - .In. N alloys and 
ordered compounds in various models: VCA (()), quadratic 
interpolation through VCA (full line), CW structures with 
Vegard’s law (A), random alloy deduced from it (dotted 
line), CW structures volume relaxed (V7), random alloy 
deduced from it (dashed line), CW with bond relaxation 
(©), random alloy deduced from it (dot-dashed line), ideal 
L1, ordering ( x ), ideal chalcopyrite (+). Ail LDA gaps are 
corrected by self-energy correction. Wurtzitic gaps are 
expected to be 0.4 eV higher. 


pairs, triplet, quadruplet. etc., correlation functions 
are considered. Here we use only the tetrahedron 
approximation. For a system in thermal equilibrium 
the P,(x,7) can be determined by minimizing the 
total free energy. However, since non-equilibrium 
growth techniques are used to deposit III-nitride 
alloys from a gas phase, it is not very relevant here 
to solve this complex problem. Since a realistic 
modelling of the non-equilibrium growth process is 
far beyond the scope of this paper, we bypass the 
problem completely by simply assuming a random 
distribution of the clusters, which implies 


P,(x) = et — xyi-* (3) 


In order to obtain the effective cluster parameters 
for the band-gap, we use the Connolly—Williams 
(CW) approach[10]. The desired quantity is calcu- 
lated from first principles for a number of ordered 
configurations (i.e. crystal structures) and is ex- 
panded in a similar manner as the disordered 
configuration in the effective cluster interaction 
parameters. If the number of ordered configurations 
considered equals the number of cluster parameters 
to be determined, the resulting set of equations can 
be inverted. This is particularly simple for the 
tetrahedron approximation within the fcc sublattice 
of cations. Corresponding to each cluster one can find 
an ordered crystal configuration in which this cluster 
is the only one that occurs. The cluster parameter can 
thus be identified directly with the corresponding 
property for that particular ordered configuration. 
Specifically, the Ga, cluster corresponds to pure 
GaN; the Ga;In cluster corresponds to the LI, 
ordered fcc based structure, in which the Ga atoms 
are placed in the face centers of the cube and the In 
atoms on the corners; and the Ga,In, cluster 
corresponds to the Ll» structure which consists of 
alternating (001) layers of Ga and In cations. The 
other ones are obviously obtained by replacing In by 
Ga. In addition to this “basis set” of CW structures 
from which we calculate the disordered random alloy 
properties, we have also considered a few other 
ordered structures to test how much the band gap 
depends on ordering. The other structures considered 
are the L1, structure which consists of alternating 
(111) layers of In and Ga, and the chalcopyrite 
(ckp) structure. The latter has been found to be 
more stable than the disordered alloy in alloys 
with large bond-length mismatch, in particular in 
In,Ga;_,As{i1]. Both correspond to 50% In. 

We first assume ideal ordered compounds, in which 
all bond-lengths are equal and assume a lattice 
constant which follows Vegard’s law, i.e. it is 
assumed to be linearly varying with the concen- 
tration. In Fig. 1, the calculated values for the band 
gap for the CW structures within this Vegard’s law 
approximation are indicated by A. The band-gap 
behavior for the random alloy deduced from it in the 
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way explained above is indicated by the dotted line. 
The chalcopyrite structure and £1, gaps are indicated 
by + and x respectively. 

We note first of all that the band-gap bowing 
within this model is slightly stronger than that in the 
VCA. Secondly, we note that the gaps of each of the 
ordered CW compounds are slightly lower than those 
of the random alloy of the same concentration. 
The gap of the (111) ordered compound (L],) Is 
remarkably lower. This indicates that layer by layer 
growth of a (111) oriented 1 + 1 superlattice may 
result in even lower band gaps for InGaN, than the 
random alloy. This is relevant because a similar result 
is expected for a (0001) wurtzitic alternating In and 
Ga layered structure. Since this is the commonly used - 
growth direction, this possibility is realistic. We also 
note that the chalcopyrite band gap Is slightly closer 
to the random alloy line than the £1) structure. This 
is consistent with the fact that the ckp 1s one of the 
so-called special quasi-random structures (SQS)[12] 
which have correlation functions close to those of the 
random configuration. In other words, this special 
structure is a good representative of the random 
alloy. 

Next, we consider relaxation effects. First, we keep 
the geometries ideal but allow the volume to relax. 
Figure 2 shows the variation of the calculated 
equilibrium lattice constant with concentration. A 
deviation from Vegard’s law is clearly visible. The 
resulting gaps of the CW structures shown by V in 
Fig. 1 and the random alloy gap deduced from it is 
shown by the dashed line. We note that this effect 
tends to increase the bandgap bowing substantially. 

Finally, we allow the individual bond lengths to 
relax. Because the InN and GaN bond lengths can 
relax towards their ideal values (those in pure InN 
and pure GaN) by a simple displacement of the N 
sublattice, only weakly coupled to the overall lattice 
constant, such a displacement is expected to occur 
fully if only bond stretch force constants were 
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Fig. 2. Lattice constant in Ga;-.In,N alloys. Vegard’s law 
(dashed line), calculated equilibrium (©) with quadratic 
interpolation (full line). 
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Fig. 3. Bond lengths in Ga; - .In, N ordered CW compounds 

calculated using Keating’s model: In-N ((]) and Ga-N(O), _ 

and random alloy bond lengths (full lines) deduced from 
them. 


present. The competition with bond angle deviations 
from perfect tetrahedral angles of 109° limits the 
bond-length relaxation. The well-known Keating 
model is adequate to deal with this optimization 
problem. The parameters for the Keating model of 
InN and GaN were derived from our calculated 
elastic constants and from experimental LO-TO 
phonon splittings[15]. From the individual bond 
lengths of the ordered compounds, the bond lengths 
in the random alloy are modelled by the cluster 
expansion 


©, Pa(x)wer* da© 
2, Pa(x)wa® + (4) 


where wS° is the number of AC-bonds in CW 
structure mn. The resulting biomodal bond-length 
distribution is shown in Fig. 3. The band gaps 
obtained for the CW structures relaxed in that 
manner are indicated by circles in Fig. 1 and the 
random alloy gap deduced from it is shown by the 
dash-dotted line. We note that this has a substantial 
effect on the band-gap bowing. 

Further work is necessary to investigate the 
combined effects of bond-length relaxation and 
volume relaxation. The bond-length relaxation could 
partially compensate for the deviation from Vegard’s 
law and thus somewhat reduce the band-gap bowing. 
Also, in our present volume relaxed approach (with 
and without cation—anion bond-length relaxation) 
the average cation—cation distance is determined by 
the equilibrium lattice constant of the corresponding 
CW structure. This means that each cation tetra- 
hedron cluster is implicitly assumed to be fully 
relaxed. Because each cluster is not independent of 
the other ones but embedded in a surrounding 
medium, this is not exactly the case in a truely 
random alloy({13,14]. We expect that neglecting this 
effect may lead to a slight overestimate of the 
band-gap bowing. 


a(x) = 


198 W.R. L. Lambrecht 


Table 1. Band-gap bowing coefficient 6 
in Gai_-In-N alloys in various models 


Model b (eV) 
VCA 0.60 
CW volume interpolated 9.89 
CW volume relaxed 1.45 
CW bond relaxed Lae 


Expt’ 1.00 


*Qsamura et ai[1]. 


The band-gap bowing is conveniently described by 
means of a bowing parameter defined as the 
coefficient of the quadratic term 6x(1 — x) in the 
alloy gap as a function of concentration. Table 1 
summarizes the band-gap bowing coefficients ob- 
tained in each of the treatments discussed above. We 
also include the experimental value for the bowing 
coefficient given by Osamura et al.[1]. 

Our calculated energies of formation £; are quite 
large, of the order of 0.3 eV/atom for the (unrelaxed) 
50% compound. The relaxation energy could not be 
calculated accurately with the Keating or ASA 
approach, but is expected to reduce £; substantially. 
Even with £; of the order of 0.1 eV/atom, though, the 
tendency towards phase separation is large in this 
system. This is consistent with experimental findings 
by Osamura ef al/.[l], who found phase separation 
behavior upon annealing above 600°C. 

In order to properly discuss energies of formation 
E, of different ordered and disordered structures, the 
chalcopyrite and Z1, structure need to be relaxed. 
This should be done by means of full-potentia! 
LMTO calculations and is postponed to future work. 
At present, for the ideal structures, we obtain 
E{L1,) < Edckp) < E{L1,). This is the expected 
ordering for lattice-matched semiconductors[!1]. 
However, relaxation effects are expected to change 
this ordering{i1]. 


4. DISCUSSION 


We find a substantially larger bowing coefficient 
than previously given in the literature[1]. However, in 
Osamura et al.[1] the concentration was obtained by 
measuring the change in lattice constant from X-ray 
diffraction data and by assuming Vegard’s law to 
hold. Because of the upward bowing of the lattice 
constant seen in Fig. 2, this procedure overestimates 
the In concentration and hence underestimates the 
band-gap bowing. Using our calculated deviation 
from Vegard’s law to correct the concentrations in 
the data of Osamura et al. we estimate that his 
measured band-gap bowing coefficient should be 
increased to ~1.8eV. We conclude that it is 
important to do an independent measurement of the 
lattice constant and the In concentration. Inhomo- 
geneities in the films, which are not unlikely given the 
tendency towards phase separation, may also 
contribute to uncertainties in the value of the band 


gap. 


Other theoretical approaches. e.g. Jenkins and 
Dow{l6] and Nakajima er al/.[17], used a virtual 
crystal approximation combined with the tight-bind- 
ing method and, as expected. obtain small bowing 
coefficients. As suggested above, some further 
refinements in the present calculations are desirable 
and may slightly overestimate our band-gap bowing. 
Nevertheless. we expect our main conclusion of a 


substantially larger bowing than previously accepted 
to hold. 

We finally note that recent reports by Singh and 
Moustakas [18] and Scholz et a/.[19} provide evidence 
of a larger bowing in In,Ga,_,.N alloys than 
previously reported{1]. Unfortunately, these reports 
do not yet allow a quantitative determination of the 
bowing coefficient. 


5. CONCLUSION 


Using first-principles LMTO calculations, we have 
shown that local variations in bonding environment, 
deviations from Vegard’s law, and bond-length 
relaxations each contribute substantially to the 
band-gap bowing parameter in Ga,_,In,N alloys, 
which is found to be about 2.6 eV. 
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’ CHARACTERIZATION OF BULK. POLYCRYSTALLINE INDIUM NITRIDE GROWN AT 
SUB-ATMOSPHERIC PRESSURES 


Jeffrey S. Dyck*, Kathleen Kash*. Kwiseon Kim’, Walter R. L Lambrecht’. Cliff C. Hayman”. 
Alberto Argoitia’, Michael T. Grossner’, Weilie L. Zhou‘ and John C. Angus? 
“Physics Dept., ° Chemical Engineering Dept., “Materials Science and Engineering Dept. 
Case Western Reserve University. Cleveland, OH 44106. 


ABSTRACT 


Polycrystalline, wurtzitic indium nitride was synthesized by saturating indium metal with 
atomic nitrogen from a microwave plasma source. Plasma synthesis avoids the high equilibrium 
pressures required when molecular nitrogen is used as the nitrogen source. Two types of growth 
were observed: 1) small amounts of indium nitride crystallized from the melt during cooling 
and 2) hexagonal platelets formed adjacent to the In metal source on the crucible sides. The 
mechanism of this latter growth is not established, but may involve transport of indium as a 
liquid film. The crystals were characterized by electron diffraction. X-ray diffraction, elemental 
analysis, scanning electron microscopy, and Raman spectroscopy. Lattice parameter and Raman 
active phonon modes are reported and compared with calculations based on the full-potential 
linear muffin-tin orbital method (FP-LMTO). 


INTRODUCTION 


Of the III-V nitride wide band gap semiconductors, indium nitride has received the least 
attention. Indium nitride is also the most difficult to grow; the equilibrium vapor pressure of N2 
over InN is higher than that over GaN and AIN [1]. Additionally, the In-N bond is the weakest 
among the [I-V nitrides (2]. High quality InN crystals are of inherent scientific interest. Also, 
for obtaining blue light emission from nitride based structures, a significant amount of In (20- 
40%) must be incorporated into GaN [3,4]. 

Efforts at bulk synthesis of InN have primarily been at high pressures because of the high 
equilibrium pressures of N2. Successful synthesis of small crystallites (5-50 micron) has been 
achieved at an N2 pressure of 18-20 kbar and a temperature range of 700-900°C [5]. We have 
recently shown that GaN and InN can be synthesized at low pressures by use of atomic nitrogen, 
N, rather than molecular nitrogen, N2, to saturate gallium or indium with nitrogen [6,7,8]. 
Subsequently, similar results were obtained by Krukowski, et al., for InN (9]. 

In this paper we report on the growth of well-faceted bulk InN using nitrogen ECR and 
ball plasma configurations at sub-atmospheric pressures. Comparisons of measured lattice 
constants and optical phonon modes are made to theoretical calculations done using the full- 
potential linear muffin-tin orbital method (FP-LMTO) [10]. 


EXPERIMENT/S YNTHESIS 


The pure indium shot was held in either an h-BN crucible for experiments without a bias 
voltage, or a machined, nitrided stainless steel crucible lined with nitrided molybdenum foil for 
experiments with bias. The reaction chamber was initially evacuated to a base pressure of 10° 
torr. In order to remove indium oxide or other impurities from the surface, the indium was 
independently heated from below while being exposed to a hydrogen plasma with 800 watts of 


microwave power at 20 torr for 30 minutes. Although the temperature at the bottom of the 
susceptor during this step was approximately 450°C. the actual temperature at the crucible was 
unknown as there is much additional heating from the plasma. After cleaning, the hydrogen flow 
was replaced with nitrogen, maintaining a pressure of 20 torr. A voltage of 20 Volts was applied 
to the crucible for bias experiments. At this point. ECR magnets were turned on for an ECR 
plasma. or left off for a ball plasma. The nitrogen pressure for ECR plasma experiments was 
lowered to 0.5 mtorr. For the growth cycle, the In was held at a susceptor temperature of 200- 
450°C while being exposed to the nitrogen plasma. 

In experiments where the molten indium metal was visible, the formation of a dull crust 
of InN on the surface of the In was observed and periodic eruptions of gas and molten metal 
occurred for higher temperatures. The gas was not identified, but was likely N>. In ball plasma 
experiments, formation of InN adjacent to the indium melt on the edges of the crucible became 
apparent after about 30 minutes. Often, most of the In was converted to InN of this form, leaving 
almost none in the crucible. After approximately 1 hour. the crucible was cooled with the plasma 
on to maintain the high nitrogen chemical potential. Additional formation of InN material 
occurred during cooling in many cases. | 


CHARACTERIZATION AND CALCULATIONS 


Scanning electron micrographs (SEM) of the polycrystalline product are shown in Figures 
| and 2. Figure | reveals growth patterns of hexagonal platelets formed on the edges of the 
crucible assembled together in “pine cone” structures. Within the clusters of platelets, the 
crystals are oriented. Typical dimensions of isolated platelets are 0.75 microns thick and 5-20 
microns wide. Additionally, InN was crystallized from the atomic nitrogen saturated indium melt 
held within the crucible. Figure 2 shows a typical growth morphology on the surface of the 
cooled indium. The “herring-bone” dendritic structure suggests uncontrolled freezing from a 
melt. 

Transmission electron microscopy (TEM) was performed on a Philips CM 20 equipped 
with Tracor Northern 5500 detector system for Energy Dispersive X-ray (EDX) analysis. A 
small fragment of material was removed from the crust formed from InN crystallization from the 
melt and glued on a Ni grid for ion milling until it became transparent for TEM observation. 
Figure 3 shows an InN crystallite. EDX analysis of the small grain with black contrast embedded 
in the InN shows a high concentration of In. We believe that it is an inclusion of In. Planar 
defects were also found in the InN grains. An electron selected area diffraction pattern was taken 
along the [1120] zone axis as shown in the inset to Figure 3. 

The indium nitride crystals were confirmed to be wurtzitic by electron and x-ray 
diffraction, and by elemental analysis. EDX analysis indicated that the only constituents were 
indium and nitrogen. In Figure 4 we show a typical powder x-ray diffraction pattern of the 
unetched material. All the major wurtzite InN lines are observed with the correct relative 
intensities. A few weak peaks due to In inclusions in the material are also seen. From the 
diffraction pattern, the lattice constants were obtained and are shown in Table I. The peak 
positions of the indium inclusions were used to correct for instrumental shift of observed InN 
peaks in calculating the lattice constants. It should be noted that our material has atomic 
spacings that are less than previously reported values. Theoretical calculation of the lattice 
parameters was done by relaxing the wurtzite InN structure by minimizing the total energy. The 
equilibrium lattice constants and elastic constants were obtained in earlier work [11]. 


Figure 1. Scanning electron microscope images of InN from growth adjacent to the indium metal source on the 
crucible sides. 


Figure 2. Scanning electron microscope image of Figure 3. Transmission electron microscope 
InN grown by crystallization from a nitrogen _ image of InN grain. Dark region is believed to be 
saturated indium melt. an In inclusion. Inset is an electron diffraction 


pattern taken along the {1 120] zone axis. 


TABLE I. Measured and calculated Lattice parameters of InN. Values are 
given In angstroms. 


| Literature Current Work 
fiz} Theorv[11] Experiment 
2) 353.36 3.525 4/-0.005 3.526 +/- 0.026 
c | 5.708...5.963 5.679 +/-0.005 ‘5.686 +/- 0.041 
cia) _1,59...1.69 1.611 +/- 0.005 1.613 


Counts (a. u.) 
Counts (a. u.) 


25 30 35 40 45 50 55 60 65 70 75 80 85 100 300 500 700 900 
2 Theta (degrees) : Wavenumber (cm-1) 
Figure 4. Powder x-ray diffraction pattern of un- Figure 5. Depolarized Raman spectrum of InN 
etched InN. Peaks are labelled with the crystallized from a nitrogen saturated In melt. 
corresponding InN reflection plane. Asterisks Narrow peaks that are not labelled arise from 
mark peaks from In metal reflections. noise at random frequencies. 


Depolarized Raman spectra were taken in backscattering geometry. Given that the InN 
material consists of randomly oriented arrays of crystallites, all six Raman active modes are 
allowed. The Raman spectrum is shown in Figure 5. The three weak peaks at 446, 475, and 488 
cm’! are likely the A; and E, splitting in the TO mode and the E) high mode. The broad peak at 
574 cm’ is interpreted as a mainly LO-type plasmon coupled mode, and the peak at 375 cm’ as 
the corresponding lower (mainly plasmon-like) mode. With this interpretation, the carrier 
density can be approximated from 


Wp = 4mne“/m’ (1) 


and using m*=0.1 to be of order 10!’ cm’. The LO-like mode shifts considerably from sample to 
sample and because of the coupling to the plasmon, the number in Table II should be considered 
an upper limit for the pure LO mode. The resulting LO-TO splitting of about 100 cm’ is smaller 
than in GaN, indicating a lower ionicity. The E: | mode was not observed. The calculated 
phonons were obtained by calculating the total energy changes and forces resulting from small 
selected atomic displacements. Assuming the harmonic approximation, these provide the 
dynamical matrix which is then diagonalized. Results of the calculations and experimental data 
are presented in Table II. The calculated (silent) B,* mode is essentially an LO-type mode at the 
zone boundary of zincblende (at L.y), folded onto the Brillouin zone center in wurtzite and thus 
should be close to the values for the LO phonons at the Gamma point. We cannot directly 


TABLE II. Phonon modes of InN. Results from the literature are also 
presented for reference. Reported values are given In cm. 


Mode Literature Current Work 
pi4j? sf 15] (16) Theory Experiment 
Ey’ 190 104 
E,° 590° 495 483 488 
ae 400 440 446 
BE 478° 490° 472 475 
aoe 596 $743 
bo 694? 
By 270 
By? 530 


* as no polarization dependence of these lines was measured, no 

distinction between A, and E, character could be discerned 

extrapolated from InGaN alloys 

“ this is more likely to be an LO mode 

4 this line coincides with a mode of sapphire which was used as 
a substrate in these experiments 


b 


calculate the latter from displacements in the unit cell because the presence of a long range 
electric field breaks the periodicity. The calculated value of B,” is consistent with our assignment 
of the experimental 574 cm’ line. It should be noted that the InN was readily converted to In2O3 
from the heating of the focused laser beam, and care was taken to avoid this contamination in the 
spectra. 


DISCUSSION AND CONCLUSION 


Bulk, polycrystalline InN is grown at pressures less than or equal to 20 torr from In metal 
and nitrogen derived from a plasma source. The method has been used to grow AIN, GaN, and 
InN without the aid of a substrate and is an attractive alternative to the high pressure techniques 
normally associated with buik growth of these materials. 

The InN has been verified to be of the wurtzite structure through both electron and x-ray 
diffraction experiments. There has been much uncertainty in the physical and optical properties 
of InN due to a lack of good bulk material to date. Existing data for the lattice parameters varies 
by about 1% and is highly dependent on the preparation technique. Our results give lattice 
constants less than previously reported in the literature and match well with calculations. In this 
study, optical phonon modes are observed through Raman spectroscopy and the fine structure of 
the TO modes is resolved for the first time. These data represent the first measurement of optical 
phonons from InN synthesized without the aid of a substrate [13-15]. Furthermore, these are the 
first ab-initio phonon calculations of wurtzite InN to the best of our knowledge. 
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